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Abstract 

The expressions for wake fields, generated in plasma (in the plasma waveguide or unlimited plasma) by 
relativistic electron bunch, was received and analyzed for the cases of the presence and absence of strong 
external longitudinal magnetic field. For the both cases the comparative analysis of the dependence of 
field amplitudes on the parameters of the electron bunch was done. 
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1 Introduction 



Presently the studies on new methods of charged particle acceleration by means of wake fields, generated 
in plasma by laser radiation (BWA (Beat- Wave Acceleration), LWFA (Laser Wake Field Acceleration)) 
and by bunches of relativistic particles (PWFA (Plasma Wake-Field Acceleration)), moving in plasma are 
intensively developed (see, e.g. reviews [1, 2] and cited there literature). The intensity of acceleration 
fields (in the order of 10 — 10 S V/ cm), attained by these methods can be used both for the charge 
acceleration, and for focusing of electron (positron) bunches in order to obtain the beams of high density 
and to ensure high luminosity in linear colliders of next generation [1, 3]. 

Linear theory of wake field generating by two- and three-dimensional rigid bunches of charged particles 
in boundless and limited plasma was developed in many works [3-11]. Nonlinear theory of wake field 
generating by a rigid one-dimensional bunch of final extent and sequence of charge particle bunch was 
developed in [12-17]. It was shown that optimum condition for wave generating is nj, = n /2 (n^, no are 
the density values of bunch and plasma electrons, correspondingly). Important result of this theory is 
the demonstration that, in the case of nonlinear wake fields a transformation ratio R = E ac /E st (E ac 
and E st are the extension of correspondingly accelerating and decelerating electric fields) depends on 
gamma-factor of accelerating bunch and may be significant without special bunch shaping, as it occurs 
in the case of linear wake fields generated by a rigid bunch. The conclusion of the results mentioned 
above is reaffirmed in [18] for the assumption (3q = vq/c = 1, (70 = (1 — /?o) = 00 ) (where vq is the 
bunch velocity), that brings to the incorrect expression for the maximum value of accelerating field E ac 
(when rib/no = 1/2 E ac = 00). 

The influence of transverse sizes of the bunch on non-linear wake field generated by short bunches 
(d <C A p , d is the bunch length, X p is the wave length) was considered in [19]. 

Non-linear theory of wake field generated by two- and three-dimensional bunches for the general case 
is not yet developed. 

This work contains the results and analysis of linear equation solution, describing the interaction 
of the axial symmetrical homogenous bunch of charged particles with plasma in the assumption of the 
plasma vorticity absence (laminar flow) , as well as at the strong external constant magnetic field applied 
along the bunch motion, when the transverse movements of plasma electrons are suppressed. 

2 Basic equations 

Vector equation, describing the excitation of nonlinear three-dimensional wake fields by the rigid bunch 
of charged particles with electron density n& moving with constant velocity Vo along the z axis through 
the cold plasma at equilibrium with density no in hydrodynamic description and in the assumption of 
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absence of the plasma vorticity 



rot (p - -^A) = 0, (1) 

is given by the following formula [3] 



i a. 



where p = p/mc is dimensionless momentum of plasma electrons, A- vector potential of electromagnetic 



field, k p — ujp/vo, uj p = ^/47rn e 2 /m-is the plasma frequency of electrons, (3 a = vq/c, and n^-arbitrary 
function of coordinates and time. 

The similar equation, describing interaction between laser pulse and the plasma, was obtained in 
[20-22]. 

Let as consider axial-symmetrical bunch, when p depends only on variable r and z = z — v a t (steady 
state). In this case vector p has only longitudinal p z and radial p r components, which do not depend on 
azimuthal angle <p, and p v = 0. 

The system of equations for the component of momentum p z and p r has the following form: 

£ (pop* vtt^) + - ^ k> : + m- (3) 

9F V J (3 av /l + p 2 - p z p n 

(3 p z -VTT7dfid. \ i d ( d Pz 



PoV 1 + P 2 -Pz dz \ rdr ) ft)V 1 + P 2 - Pz rdT V dT 

Pz 19/9 



1 9 2 /9r 



^(^1+?-*) + ^ (4) 



9,z \ r 9r 



r 9r ( 9r ^ ^ ) } 



Inserting in (3) and (4) p r = and p 2 = /O z (z) (dependence on r is absent) we come to the one- 
dimensional nonlinear equation, considered in [12-17] 

^ (PoPz - VTTpI) + - f *^ + «^ = 0. (5) 
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Inserting in (4) p r = and expressing the derivative dp z /dz through derivatives of p z we receive for 
scalar potential ip the following equation, describing interaction of electron bunch with plasma in the 
presence of the external magnetic field B (0, 0, Bo) [19]: 

1 d ( d X \ , 1 d 2 X , fl2;2 (, PoX \ _ Potf n b 



r- 



^4 + ^1- , == 5 =^/-, (6) 



rdr\dr) ' 7o 2 ^ ' " u '" p \^ Vx 2 - lR J 1% n o' 

where 



mc 2 7q 



v/l + pf-A/V (7) 



3 Wakefield generating at vorticity absence 

Let us consider the problem of wake field generating by rigid cylindrical bunch of radius a and horizontal 
dimension d with homogeneous distribution of electrons of the bunch 

0<r<a, 0<z<d, 
n b ={ _ " Z ( 8 ) 

a <r <b, z > d, z < 0, 

moving in conducting plasma waveguide of radius b> a. 

Assuming in (3) and (4) n b /no <C 1 (linear approximation) and linearizing the system on p we shall 

obtain the following system of equations, describing the process under consideration: 

dzdr + rdz dr 2 r dr + ^ dz* + ^> Pz ^>n ' (9) 

a% 2 2 _ 

ffzdr 7o 2 dz 2 + ^> pr - U ' 

To define p z (r, z) and p r (r, z) let us perform the Hankel transformation [23] of the equations system (9) 
on r in the finite limits (0, b) and solving the obtained equations on z in the assumption of the continuity 
conditions of the momentum components p z , p r and components of the electrical field E z = mc e V(> ^fc , 
E r = mcv ° ^te at the front (z — d) and rear (z = 0) bunch boundaries we shall receive the following 

e dz 

expression for the field components E z and E r : 
Inside the bunch (0 < r < a, < z < d) 



El = 



mv uj p n b 



n 



sin k p (d — T) x 



(10) 



x < 1 - (k p a) 



Io(k p r) 
Io(k p b) 



[K (k p b)Ii(k p a) + ^(kp^K^kpa] 
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j mv uj p n b ~Ji(k p r) 

E r = (k p a) cos k p (d ~ z) 1 x (11) 

e n h{Kpb) 

x [h(k p a)K (k p b) + IoikpfyK^kpCi)} + Ej. H ; 



Over the bunch (a < r < 6, < z < d) 



mv Lo p n b . 7i(fc p a) 



e n lo(kpb) 
x [K {k p r)I {k p b) - I {k p r)K (k p b)} + E 1 ^, 



mvoUjpn b h (k p a) 



E 1 / = ^ — (fc p a)cosfc p (d-?) ; ,\ x (13) 

e no lo{kpb) 



x 



where 



^ - 2 ^^^T g « (ip + ,1) JIM &h70K U - z ) ' (14) 

in the corresponding areas on r. 

Before the bunch (0 < r < a < b, d < z < oo) wc have the following expression for the field 
components: 



HI _ mvoujp nb (k p a) ^ M „J (gg£) ^ ( / d\ 



(16) 



After the bunch (— oo < z < 0) we have 



El y = !jm^ L !}± [sin{kpI) + sinkp{d _ I)]x (lg) 

e n 



|l - M^^j [^o(fcp6)/i(Apa) + Jo(M)^i(M)] } 



I jpIV 



r<a, 
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e j V = mtwpn, [sin(fcp?)+sinfcp(d _ ?)] x (lg) 
e no 

x(fc p a)^TTT [VMWVO - K (k p b)I (k p r)} + E%, 
a < r <b, 



where 



(20) 



*T = -2^=-^(fc pa )sm[ ^\ S mk p \z-^) x (21) 



? iv o mv u p n h / fc p <A . , (~ d 



x 



[/i (fc p a)X (k p b) + I (kpb)^ (k p a)} , atO < r < a 



[K (k p b)h(k p r) + 7 (fc P 6)^i(fc P r)] , ata < r < b 

d 



The components different from zero of _E-wave (-B z , £7 r , _B V ) magnetic field is determining by the 
formula 

mc 2 ( dp r dp 

and is described by the following expressions 



B * = — (22) 



nl II r ,mv uj p ,n b Ji (pnl) Ji (p n f) 

B P = - 



n ° K 2 b 2 J{(p n ) 



! - i 70 2 cIi7qk ( — — z 



x 

atO < 2" < d, < r < 6, 



B - - -2^/?o(fc P «)-^ K 2 52Ji2(Mn) e U J&h ^°Tj' (24) 



d < z < oo, < r < 6, 



n=l 

— oo < z < 0, < r < 6. 
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In formulas (10)-(25) Jo, Ji are the Bessel functions, and I , I\, K , i^i-are the modified Bessel 
functions, k 2 = k 2 (3 2 + [i 2 Jb 2 . 

As it may be seen from the given expressions for field components E z , E r they consist of periodic 
(wake field) and non-periodic ("Coulomb") parts. Before the bunch the field has only non-periodic part, 
which exponentially falls of with the remote from the front boundary, so that part may be neglected. After 
the bunch the " Coulomb" part also falls off exponentially with the remote from the rear boundary z = 
and only periodic wake field remains. Magnetic field component B v has only non-periodic "Coulomb" 
part, which exponentially also decreases at the remote from the bunch boundaries and it can be ignored. 

In the range of < z < d, where magnetic field is not small, the radial force f r = —eE r + eft^B^ 
acting on the bunch electrons, and this force in some ranges on z can compress the bunch, focusing it. 
Before the bunch (z > d) and in the range of wake field (z < 0) the component of the magnetic field B v 
is small and the radial force is f r —eE r . 

Below one can find an expression for the radial force f r acting on the bunch (0 < z < d) at (3 ~ 1, 
7o > 1: 



-mWpUQ — (k p a) [1 — cos k p (d — z)] x (26) 

^^[hik^Koik^ + Ioik^K^a)}, r<a, + ( ±_ 

[K (kpb)h(k p r) + Io(fcp6)tfi(fcpr)] - a<r<b, W 

where O^q 2 ) is the smaller defocusing part of the force. 

Let us also give an expressions for the fields E z and E r inside and after the bunch in the case when 
a = b and in the case when b — > oo (unlimited plasma) : 



K 

El 

E 1 / 
a 



mvQLOp ni, 
e n 
mv LUp 71b 
e n Q 



sin k p (d—z) 
cos k p (d — 7) 



_ Io(k p r) 
I (k p a) 
h{k p r) 
I Q (k p a) 



E: 



h- 



— [ s in k p (d — z) + sin (k p z)] 
[cos kp(d — Z) — cos (k p z)] 



e n 
mvouip rib 
e n Q 



1 - 



Io(k p r) 
I (k p a) 
h(k p r) 
I (k p a) 



+ E IV 



E 



IV 



b, 



(27) 



where 



E - - 2 — ~, To g ^MWTJ^ e ' Bbwi U " ') • <28) 
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j _ mvoujp n b ^ Ji ( Mk £) 



70 2"ch7 K ( ^ 



and ^ 



at z < 0. 

At the a — ► oo I (kpCi) — > oo and the expressions for £^ and £^ y coincide with the expressions for 
one-dimensional bunch, while E\. and E^ v — > 0. At fc p a <C 1 (fc p r < k p a <C 1) periodic parts of the fields 
El ,IV and El' IV are proportional to (k p a) 2 and fc p r. Non-periodic parts are small at 70 ~> 1. 

The expressions for £7 2 and £7 r inside the bunch and in the wake field at b — > 00, 70 3> 1 go over into 
corresponding expressions for the case with unlimited plasma. In this case at k p r <C 1, k p a <C 1 (r > a 
or r < a) the longitudinal fields inside and over the bunch as well as in the wake are proportional to 
(k p a) 2 , while the radial components of the fields E^ IV <~ r/2a at r < a and El' IV ~ a/2r at r > a, 
e.g. they increase with the remote from the bunch centre inside the bunch and decrease with the remote 
from the bunch boundaries inside it. 

Let us define the plasma density for the mentioned four regions. From the Poisson equation it follows 
that 



n K 



1 

Aire 



8E E r 8E Z . 
+ — + - 4-7re (n - n b ) 



dr ' r ' dz 

In the regions outside the bunch in (29) we should suppose n& = 0. Using the expressions (10)-(21) for 
the fields we shall find out the following expressions for the plasma densities: 



n 1 — ^ (1 — cos k p (d — z)) , r < a, 



< z < ± 



no, 



ni 11 = n , 



a < r < b, 
atO < r < b, d <z < 00, 



(30) 
(31) 



n 



1 



no 



(cos (kpl) — cos k p (d—z)) 



00 < z < 0. 



(32) 



r < a, 

no, a < r < b, 

Thus, the density n e depends periodically on z inside and after the electron bunch at r < a and does 
not depend on r and sizes a and b. 

For the regions outside the bunch a < r < b and < z < 00, —00 < z < n I J^ IV = n and 
the density n e undergoes change to the lateral surface of the bunch and in the wake at r = a. On the 
waveguide surface r = b, n e = uq. 



4 Wake field excitation under strong external magnetic field 

The linear equation for the potential <p, describing the interaction of cylindrical electron bunch of radius 
a and length d with unlimited cold plasma, follows from expressions (6), (7) in assumption ef/mc 2 <C 1, 
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nb/no <C 1 and has the following form: 

a/ 2 ? 1 8^ t i aP<p , fcp 



2 2 mn h 



(33) 



dr 2 r dr ^ dz 2 7 2 
where n& is given by the expression (8). 

After making Hankel (Fourier-Bessel) transformation of the equation (33) on r in boundless limits 
(0, oo) under the condition that tp(z,r — > oo) = we shall receive the following equation 



d 2 ip{a,z) 2 ^ 



dz 2 



— A ip(a, z) = h, 



(34) 



where f(a,z) = f °° <p(r,z)J (ar)rdr, A 2 = j 2 (a 2 - fc 2 ), k 2 = k 2 /^ 2 , h = ^k 2 v 2 ^-f 2 ^Ji(aa), J and 
Ji-are the Bessel functions. 

Before the bunch (z > d), where n b = 0, h = the potential tp(r, z) is assumed as equal to zero (we 
ignore the "Coulomb" field, see chapter 3). 

Solving the equation under assumption of continuity of the potential Tp{a, z) on the front (z = d) and 
rear (z = 0) boundaries of the bunch, we shall come to the following expressions for ip(r, z) in the ranges 
inside and over the bunch (0 < z < d, < r < a, a < r < oo) 



<Pi 



— kzvf. — aRe 
e p n jq 



Jo(ar) Ji (aa) 



1 - e 



-\(d-z) 



da, 



a 2 — k 2 

where k has the positive imaginary part Imfc = fc > 0, A = 7oVa 2 — k 2 at a > k, A = 
a < k. 

After the bunch (z < 0, < r < oo) the potential y> has the following form: 



(35) 

—ijoVk 2 — a 2 at 



m 2 2 -.,. ,-. 

(,52 = fc^O a ^ C 

e ' 



r Jo (a 
io a 



Jo(ar)Ji(aa) 



da. 



(36) 



' n ./n a 2 - k 2 

Because the " Coulomb" components (a 2 > fc 2 ) are small one can ignore them and bring out the expression 
of square brackets from the integral in the point a = 0, where it makes the main input into the integral. 
In this case the expressions (35), (36) significantly simplify: 



m 2 2 n b 

(pi = k v — aRe 

e 1 n 



1-e 



ik p (d—z) 



Jo (or) Ji (aa) 
a 2 — fc 2 



da, 



<f 2 = k n v n — aRe 

e p n 



-ikp(z-d) 



Jo(ar)Ji(aa) 



f°° J (ar)J 
Jo a 2 - 



k 2 



da. 



The components of the fields are determined from the expressions: 



E x = - 



1 dip 



7o dz 



.— , E r = — 



d<p 
dr ' 



Bin — PoE r . 



(37) 
(38) 

(39) 



Inside the bunch (0 < r < a, < z < d) we have: 
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El(r,z) = 1 ^^—\smk p (d-zl\l + l(ka)J Q (kr)Y 1 <ka) 
e no L 2 

— — (ka) cos k p (d — Z) J (kr)Ji(ka) j , 



(40) 



niVQLOp 1 Tlf, f IT 



e 70 n 



|| [i _ cosfcp (d - I)] (jfeo)Ji(fcr)yi(fea)- (41) 



~(fca) sinfcp (d — z) Ji(fcr) Ji(fca) j , 



where fc = fc p /7o- 

Over the bunch (a < r < oo, < z < d) the field components are given the following expressions: 



mvoUp rib ( it 



e n 



j|(fca) sinfcp (d — z) Ji(ka)Y (kr) — 



-|(fca) cos k p (d—z) J (fcr) Ji (fca)| 



El^r.z) = H^P— ^^(ka) [1- cos k p (d- I)} J 1 {ka)Y 1 {kr)- (43) 



e 7o "-o 



— |(fca) sin kp(d-z) Ji(fcr) Ji(fca) j . 
The wake field is determined by the following expressions: 



Ei v (r,z) = 



mvQLOp rib 
e n 
+ [cos 



< r<a, 



mvQUjp rib 
e n 



| [sin (fcpz) — sinfcp (z — d)] —(ka)Ji(ka)Y (kr)+ 



7T _ "I 

+ — (fca) [cos (k p z) cos hp (z- d)],h{ka)J Q {kr)\ 

a < r < oo, 



(42) 



| [sin (fc p 2) - sinfcp (z - d)} 1 + ^(ka)J (kr)Y 1 (ka) + (44) 
(k p z) — cosfcp (z — d)] —(ka)Jo(kr)Ji(ka)^ , 



(45) 



El v {r^ = ^^-l^|[cos(fcp?)-cosfcp(2 r -d)]|(fca)J 1 (fcr)Y 1 (fc a )- 
— [sin (fcp?) — sinfcp (z — d)] — (ka) Ji(fcr) Ji(fca) j , 



r < a, 



(46) 
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E™(r,z) = V±{ [cos (k p z)- cos k p {z-d)} J(Ao)Ji (fca)Yi(fcr)- (47) 

e 70 no L 2 

— [sin (k p z) — sinfcp (z — d)\ —(ka)Ji(kr)J\(ka) j , 
a < r < oo. 

The components of the magnetic field B v is determined from the expression B v = (3oE r . 

The comparison of the expressions for the components of the fields (10)-(25), generated by the electron 
bunch in plasma without magnetic field, with the corresponding expressions (40)-(47) with the constant 
longitudinal strong magnetic field B shows, that in the last case the character of the fields qualitatively 
varies from the case of plasma without the field. On the one hand, the dependence on the transversal 
coordinates r and a is determined by the Bessel functions Y and J, having the oscillating character, 
and on the other hand, there is more distinctly expressed dependence on 7-factor of the bunch. Besides, 
magnetic wake field B v is equal to zero in plasma without the field and defers from zero in the case of 
plasma with B 7^ 0. 
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